Parton distributions and small-x QCD at the Large 
Hadron Electron Collider 

Juan Rojo and Fabrizio Caola 

Dipartimento di Fisica, Universita di Milano and 
INFN, Sezione di Milano, 
Via Celoria 16, 1-20133 Milano, Italy 

The proposed Large Hadron Electron Collider (LHeC) at CERN would bring Deep- 
Inelastic scattering into the unexplored TeV regime. The LHeC rich physics program, 
among other topics, includes both precision SM measurements to complement LHC 
physics as well as studies of QCD in the high energy limit. The present contribution 
reports on ongoing studies within the NNPDF framework towards the LHeC CDR. We 
study the impact of LHeC simulated data on PDF uncertainties, in particular the small- 
X gluon. We also assess the LHeC potential to disentangle between various scenarios 
of small-x QCD, including saturation models and small- a; resummation. Finally, we 
explore how deviations from DGLAP can be quantified in inclusive measurements. 

Introduction The Large Hadron Electron Collider (LHeC) [T] is a proposal for a Deep- 
Inelastic scattering facility in the TeV range which would operate in parallel with the LHC. 
It would use the 7 TeV LHC proton beam colliding with a high energy electron beam, coming 
either from a LHC-like ring or from a linear accelerator. The kinematical coverage of such 
machine would extend the HERA kinematical coverage by two orders of magnitude both in 
X and in Q^. 

After the experience at HERA, it is clear that the LHeC physics potential includes 
the capability to probe the nuclcon structure and its flavour decomposition with very high 
precision. However, the standard approach to PDF determination [51 [3] suffers from sev- 
eral shortcomings. The most important ones arc related to the fine-tuning of the PDF 
parametrizations and the statistical definition of the associated PDF uncertainties to the 
available dataset. 

These shortcomings render difficult its application to extrapolation regions like the LHeC 
kinematics. In particular, in the standard PDF approach, PDF uncertainties are artificially 
reduced in extrapolation regions due to relatively simple polynomial parametrizations em- 
ployed, thus making difficult a quantitative assessment of the impact of new data from 
unexplored regions into the PDFs. On top of that, subtle deviations from DGLAP evolu- 
tion which might be present at small- a: are difficult to probe with simple fixed functional 
forms because their lack of flexibility could lead to misleading results. 

A method to bypass the above problems has been proposed by the NNPDF collaboration. 
Within the NNPDF approach [H [H El El [HI (3 (see also [TO])! a combination of neural networks 
as universal unbiased interpolants with Monte Carlo sampling of experimental data for error 
propagation render the PDFs and associated uncertainties statistically faithful0 

In this contribution wc report on ongoing studies of PDF determination and small- a; 
QCD within the NNPDF approach. In particular, we concentrate on LHeC pseudo-data 
at small-x. We consider F2{x,Q^) and Fl{x,Q^) simulated pseudo-data at small-x, in a 
scenario in which the LHeC machine has electron energy of Ep = 70 GeV and electron 

The NNPDF methodology has also been applied to other physical problems in 1121 [13] . 
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Figure 1: The kinematical coverage of the LHeC pseudo-data used in the present studies, together 
with the data already included in the reference NNPDF1.2 dataset. 



acceptance of 9^ < 179°, for an integrated luminosity of J C = 1 fb~^. Full NC and CC 
cross-section simulated pseudo-data in various other machine scenarios are available, and 
their are under current scrutiny. 

The reference baseline for the studies presented in this contribution is the recent NNPDF1.2 
parton set [S] , a PDF analysis of all relevant inclusive DIS data together with neutrino charm 
production to constrain strangeness. The kinematics of the pseudo-data, together with that 
of the NNPDF1.2 analysis are shown in Fig. [1] The average total uncertainty of the simu- 
lated F2 pseudo-data is ~ 2%, while that of Fl is ~ 8%. 



Constraining parton distributions at small-x In spite of the wealth of precision data on 
small- a; structure functions at HERA, some PDFs, most notably the gluon, have still rather 
large uncertainties in this region [7]. In order to quantify how these PDF uncertainties 
would be reduced with LHeC data, we have repeated the NNPDF1.2 analysis with the 
addition of the LHeC pseudo-data, with central values from the NNPDFl.O predictions and 
experimental uncertainties corresponding to the simulated LHeC scenario described above. 
The joint data set is shown in Fig. [1] 

First of all, we include only F2 LHeC pseudo-data into the analysis. Although the fit is 
as expected perfect, the reduction of the small- a; uncertainties is rather moderate, as shown 
in Fig. [2] (left). Our results therefore indicate that only F2 data, even if very accurate, is not 
enough to pin down the gluon at small-x, due to the fact that the gluon PDF only enters 
through scaling violations and higher order corrections. 

The next step consists of the addition of the complete F2 and Fl pseudo-data. In this 
case, the joint fit with F2 and Fl pseudo-data leads to a sizable decrease of the small- a; 
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Figure 2: Left: the gluon -which results from a fit -which includes on top of the NNPDF1.2 dataset 
[9] only F2 LHeC pseudo-data. Right: the gluon PDF determined from the same dataset together 
-with F2 and Fl LHeC pseudo-data, generated from the central and ±l-o" NNPDFl.O gluons. 

gluon uncertainties, -which can be understood from the greater sensitivity of Fl to the gluon 
PDF. To quantify more these results, "wc have generated LHeC pseudo-data in three different 
scenarios: one "where the gluon is the central NNPDFl.O gluon and t"wo more "were the gluon 
sits near the associated ±1-ct envelope. We have repeated the joint F2 +Fl analysis in these 
three cases: results are sho"wn in Fig. [5] (right). As expected, after the fit the three extreme 
scenarios for the small-x gluon can be precisely disentangled. Therefore, it is clear that the 
LHeC has the potential to pin do"wn "with great precision the behaviour of the lo"w-a; gluon, 
but only after accurate measurements of Fl are performed. 

Probing small-x QCD The LHeC "would also provide us "with an improvement of our un- 
derstanding of the small-x dynamics of QCD. Indeed, even after years of intensive study at 
HERA, no convincing evidences for departures from standard DGLAP evolution has been 
found. For example, geometric scaling of HERA data, "which "was thought to provide a 
clear signal for saturation, "was recently shown to be consistent -with linear QCD evolution 
as -well [14]. The situation could be different at the LHeC, -with its extended kinematical 
coverage at small-x (see Fig. [1]). 

In order to test -whether or not a DGLAP analysis can reproduce theoretical predictions 
-which deviate from pure DGLAP in inclusive measurements, LHeC pseudo-data has been 
generated not -within the DGLAP frame-work, as in the previous sections, but rather from 
t-wo different models: the AAMS09 model [T^, -which is based on BK evolution -with running 
coupling, and the FS04 model [TS|, based on the dipole model. 

We have repeated the PDF analysis of the previous section but -with these ne-w pseudo- 
data. Although clearly the procedure is not consistent (for example, PDF error reduction 
-would be meaningless in this case) , it provides an illustration of a potential analysis technique 
-which ultimately should be applied to experimental data. For both the AAMS09 and the 
FS04 models the conclusions of the study are the same: the DGLAP analysis reproduces 
perfectly the F2{x,Q^) pseudo-data, -which implies that although the underlying physical 
theories are different, from a practical point of vie-w the small-x extrapolations of AAMS09 
and FS04 for F2 are rather similar to DGLAP-based extrapolations. 
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The situation however is different for Fl{x,Q^): provided the level arm in is large 
enough, the DGLAP analysis fails to reproduce simultaneously Fl in all the bins, and 
thus the overall is very large, a clear signal of the departure from fixed order DGLAP 
of the simulated pseudo-data. This effect is illustrated in Fig. [3l where the results of the 
DGLAP analysis are compared with the LHeC pseudo-data generated from the AAMS09 
model. 

There exists however other scenarios for QCD at small- a; than saturation/dipolc mod- 
els. In particular, linear QCD evolution with rcsmimration of BFKL small- .t logarithms 
is the natural extension of standard DGLAP evolution. Recently, the full set of small- a; 
resummed splitting functions and coefficient functions became available |17| . The results 
of Ref. [TT] have been used to compute resummed K-factors [TB], defined as the ratio of 
structure functions NLO small- a; resummed over fixed order NLO, as a function of (a;, Q^). 
These K-factors can be used for realistic, though qualitative, phenomenological studies of 
the impact of small- a; resummation. 

We have used these resummed K-factors to estimate the feasibility of the LHcC to dis- 
entangle between scenarios for small-a; linear QCD: NLO, NNLO and NLOres. In Fig. [4] we 
show the LHeC pseudo data for F2{x, Q"^) at small-x compared with the NLO NNPDFl.O 
prediction (including the associated PDF uncertainties) and the corresponding NNLO and 
NLOrcs computations, obtained from the NLO one with these K-factors. Fig. [4] seems to 
indicate that a PDF analysis capable of implementing both the the NNLO and NLOres 
computations of physical obscrvablcs has the potential to disentangle between these two 
scenarios of small-a; QCD, given the foreseen experimental accuracy at the LHeC. 

Departures from DGLAP It is clear from the previous discussion that there is some con- 
tradiction between the two goals of our study: either we determine the PDFs or we find 
evidence for saturation or resummation. However, both these goals require the same first 
step: we have to determine the kinematic region where saturation/rcsummation effects, or 
more general, departures from fixed-order DGLAP evolution, start to play a role, if any. 

The idea is therefore to single out a safe region, where the standard PDFs extraction via 
fixed order DGLAP is reliable, and a small-x region where deviations from pure fixed order 
DGLAP could provide evidence for saturation or resummation. The determination of these 
kinematic regions is a highly non trivial task: both BFKL and non-linear effects are known 
to be rather moderate in the HERA region, and thus are difficult to observe in inclusive 
observables. In particular, they could be absorbed in the initial condition for ffexible enough 
parametrizations of the PDFs. This might already be the case at HERA for i^2, and if so 
even more at the LHeC. 

A possible approach to this problem is the following. First we repeat the global PDF 
analysis removing subsets of data where small-x effects could play some role. Then we de- 
termine whether NLO DGLAP is able to reproduce the excluded data or not. A tension 
between the actual data and the DGLAP prediction should mark the onset of some satura- 
tion/rcsummation effect. As a cross check, we can assess the NLO DGLAP fit equality in the 
fitted data region: the fit quality should improve if there is some tension between DGLAP 
and the actual data in the excluded region. Note that deviations from DGLAP are known 
to be rather moderate, hence our approach is meaningful only on statistical grounds. It 
is therefore mandatory to perform a PDF analysis with no parametrization bias and with 
faithful uncertainty estimation [7j . Related studies of the stability of global analysis within 
the standard PDF approach have been reported in [19l [20] . 
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Figure 3: The results of the combined DGLAP analysis of the NNPDF1.2 data set and the LHeC 
pseudo-data for Fl{x,Q^) in various bins generated with the AAMS09 model. 
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Figure 4: A comparison of various approximations to linear low- a; QCD for F2 at the LHeC: the 
NNPDFl.O prediction which includes PDF uncertainties (green lines) and the NNPDFl.O result 
corrected with the NNLO (black, dot-dashed) and NLOres (violet, short-dashed) K-factors. The 
expected experimental precision at the LHeC is also shown for illustration. 
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Figure 5: Left: distances, Eq. [T] in the < 1.5 x""'^ HERA region computed using a global fit 
with these data points included. Right: distances in the same region computed when these points 
are excluded from the global fit. 



This approach has been apphed to search for DGLAP deviations in the small- a; HERA 
data. Taking as a reference the NNPDF1.2 analysis [S], we excluded data points with a 
saturation-inspired cut > Ql{x) = Ax~^-^, with A ranging from 0.2 to 1.5. In order to 
quantify deviations from DGLAP we computed the distance d{x,Q^) between the DGLAP 
extrapolation F^^ for an observable F and the actual data F'^^*^'^, defined as 

V '^Rt + '''data 

A typical result for the cut > 1.5 x"" '^ is shown in Fig. [5l Note that while in the 
global fit distances seem uncorrelated, in the fit with the kincmatical cut > 1.5 a;~° '^ 
there seems to be a hint of a correlation, that is, the NLO DGLAP prediction tends to be 
smaller than actual data. A systematic study is in progress using these methods in order to 
determine the statistical significance, if any, of departures from DGLAP in inclusive small- a; 
data. The ultimate validation of the method will be its application to LHeC pseudo-data, 
where there the underlying physics can be varied within various scenarios. 



(1) 



Outlook This contribution summarizes some of the studies performed within the NNPDF 
framework in order to assess the physics potential of the LHeC as a probe of the nucleon 
structure and of small- a; QCD dynamics. From these preliminary studies, one solid conclu- 
sion is that the importance of accurate measurements of Fl(x,Q^) should be emphasized. 
Ongoing work towards the LHeC Conceptual Design Report includes the generalization of 
the PDF analysis to the complete LHeC data set for various scenarios and the impact of the 
reduction in PDF uncertainties on LHC phenomenology. 
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